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' Abstract. Production of final states containing open charm (c) and beauty (b) quarks at HERA 

l/^ I is reviewed. Photoproduction (PHP) of the charm meson resonances D*, and Ds, as well as 

D* production in the deep inelastic scattering (DIS) regime, are measured and compared to QCD 
predictions. The excited charm mesons i:'5'(2420), (2460) and £ifi(2536) have been observed 
and the rates of charm quarks hadronising to these mesons were extracted. A search for radially 
excited charm mesons has been performed. PHP and DIS beauty cross sections are higher than 
1^ ^ expected in next-to-leading order (NLO) QCD. 

o 

INTRODUCTION 

^ . The HERA e-p collider accelerates electrons (or positrons) and protons to energies 

(D : of Ee = 27.5 GeV and Ep = 920 GeV (820 GeV until 1997), respectively. The HI and 

Q_i' ZEUS experiments are located at two collision points along the circulating beams. The 

^ ■ incoming interacts with the proton by first radiating a virtual photon. The photon is 

either quasi-real with < 1 GeV^ and Q^^dmn ~ 3 ■ 10"^ GeV^ (PHP regime) or highly 
virtual (Q^ > 1 GeV^ - DIS regime). 
^ ■ The large masses of the heavy quarks (HQ) c and b provide a "hard" scale needed 

for the comparison of data to QCD predictions. In leading-order (LO) QCD, two types 
of processes are responsible for the PHP of HQ's: Direct photon processes, where the 
photon interacts as a point-like particle with a parton from the incoming proton, and 
resolved photon processes, where a parton from the photon partonic structure scatters 
off a parton from the proton. Heavy quarks (Q) present in the parton distributions of the 
photon lead to LO resolved processes, such as Qg — > Qg (where g is a gluon), which are 
called heavy flavour excitation. In NLO calculations, only the sum of direct and resolved 
processes is unambiguously defined. 

Two different NLO calculations are available for comparison with measurements of 
HQ PHP at HERA: 1) a fixed-order ("massive") approach [1, 2], where HQ's are pro- 
duced only dynamically in the hard subprocess. This calculation is expected to become 
less accurate when pj_ >> niq, where p±_ and niq are the transverse momentum and 
mass of the HQ; 2) a resummed ("massless") approach [3, 4], where the massless HQ's 
from the photon and proton parton distributions are used explicitly. This calculation is 
expected to yield better results as >> mg. 
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FIGURE 1. (a) M{KTnTs) — M (Ktt) distribution in the _D" mass region (full dots). The histogram is 
the distribution for wrong charge combinations, (b) M{K^tt^) distribution for events excluding the D*^ 
region. The solid curve is a fit to two Gaussian shapes for the right and wrong K mass assignments plus a 
sum of exponential and linear backgrounds, (c) M{K^K^n^) distribution for events inside the mass 
range (1.0115 < M{K'^ K~) < 1.0275GeV). The solid curve is a fit to a Gaussian resonance plus an 
exponential background. 



PRODUCTION OF D ^, AND Ds MESONS 

The charmed meson D*"^ has been reconstructed via its decay chain D*^ 
D^Tig — > {K^n^)TXg (+ C.C.). Fig. 1(a) shows the mass difference distribution, 
AM = M{K-K-Ks) - M{K7t), in the D° mass region 1.83 < M{K7t) < 1.90GeV 
in the kinematic range p^* > 2 GeV and \r]^*\ < 1.5, where p± is the transverse mo- 
mentum and 1] = — lntan(6'/2) is the pseudorapidity. The polar angle, 9, is defined 
with respect to the proton beam direction. The plot includes PHP and DIS ZEUS data 
collected during 1995-2000 [5]. A clear D*"^ signal is seen (dots) on top of a small 
combinatorial background, estimated by wrong charge combinations (histogram), where 
both tracks have the same charge and its has the opposite charge. Defining D*"^ 
candidates as events with 0.144 < AM < 0.147 GeV, a signal of 31350 ±240 D*^ 
mesons was found after background subtraction. 

Inclusive production of charm hadrons other than the D*"^ have also been observed. 
The M{K^n^) distribution for PHP events excluding the D*^ region (0.143 < AM < 
0.148 GeV) is shown in Fig. 1(b) for a restricted kinematic region, using a ZEUS 
event sample with integrated luminosity C = 66 pb^^ [6]. A clear signal is seen. 
The production ratio, Py, of vector to pseudoscalar-i- vector ground state (orbital an- 
gular momentum L = of the cq system) charm mesons can be approximated by 
= {a{D°)/a{D*+) - BD*+^DOny\ where a{D°) and a{D*+) are, respectively, the 
inclusive and D*"^ cross sections and Bd^+^^o^ is the D*^ branching ratio to D^n^. 
Using also the D° signal from events in the D* region with the same cuts yields the pre- 
liminary result Py = 0.546 ± 0.045(stat.) ± 0.028(syst.). This measurement is in good 
agreement with the e+e~ annihilation results, 0.57 ±0.05 and 0.595 ±0.045 [7], sup- 
porting the universality of charm fragmentation. 

Using a ZEUS PHP sample with 38pb~^ in a restricted kinematic region, a Dg 
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FIGURE 2. (a) Differential D* cross section in logio(Q^) compared to the NLO HVQDIS calculations 
with different structure function parametrisations. (b) Differential D* DIS (1 < < 100 GeV^) cross 
section in 77^ compared with HVQDIS (lower shaded band) and CASCADE (upper shaded band) 
predictions, (c) Same as b) for 4 < < lOGeV. 



signal is seen (Fig. 1(c)) in the K+fT^vr^ mass distribution for events where 
M{K'^ K~) is in the (p mass range [8]. The ratio of the Dg to D* cross sec- 
tions in identical kinematic regions was measured to be cTep^DsX /o'ep^D*x = 
0.41 ± 0.07(stat.)+o;o5(syst.) ± 0.10(br.), where the last error is due to the uncer- 
tainty in the Dg (pn branching ratio. This result is in good agreement with the ratio 
/(c^ D+) / f{c^ D*+) = 0.43± 0.04± 0.11 (br.) obtained [8, 9] from e+e" 
experiments, again confirming the universality of charm fragmentation. 

D*^ Production in DIS 

Open charm production in the DIS regime is dominated by boson-gluon fusion (BGF) 
processes, where the boson (gluon) is emitted from the incoming electron (proton). 
Fixed-order NLO perturbative QCD (pQCD) calculations are available in the form of 
a Monte Carlo (MC) integrator (HVQDIS) [10]. The ZEUS preliminary D* differential 
cross section in Q^, using a sample of £ = 82.6pb^^, is shown in Fig. 2(a) for the 
kinematic region Q"^ > lOGeV^, pf > 1.5 GeV and |?7^'| < 1.5 [11]. The distribution 
compares well with the HVQDIS calculations, using 3 different parton distribution 
functions in the proton and a charm quark mass range 1.3 < rric < 1.6GeV. For Q"^ >> 
ml, resummed NLO calculations should be superior. However, up to ~ 1000 GeV^, 
the data is nicely described by the fixed-order scheme. The D*^ cross sections were also 
measured separately for e+ and beams [12]. Integrated over > 20GeV^, the e^p 
cross section is higher than that for e^p by ^ 3 standard deviations. Both results are 
compatible with the predictions within the theoretical uncertainties. 

The HI Collaboration has measured D*^ production in DIS, using C = 18.6pb^^, 
in a kinematic region 1 < < lOOGeV^ p^' > 1.5 GeV and jr^^'l < 1.5 [13]. The 
differential cross section in r]^* is compared in Figs. 2(b-c) to HVQDIS [10] and 
to CASCADE [14] calculations, which implement a version of the CCFM evolution 
scheme [15]. The shaded bands reflect the uncertainties in the predictions due to rric 
(1.3 — 1.5 GeV) and the allowed fragmentation parameters. CASCADE shows better 
agreement with the overall r]^* distribution (Fig. 2(b)), while the HVQDIS prediction is 
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FIGURE 3. Differential cross sections in and 77^, where D stands for D* or Dg- The Dg (dots) and 
£)* (triangles) data are compared with NLO predictions for Dg (full curves) and D* (dashed curves) with 
two parameter settings: = l.SGeV, /iji = m± (thick curves) and = 1.2GeV, jiji — 0.5m± (thin 



curves). Here is the renormalisation scale and mj_ — \fm^ 
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too low for the forward region. CASCADE agrees poorly with the data at the low-rj 
high-pf region (Fig. 2(c)). 

Photoproduction of charm mesons 

Differential Df cross sections in p± and r] [8] were compared to those for D*"^ [16] 
in the same kinematic region < 1 GeV^ 130 < W < 280 GeV, > 3 GeV and 
|?7^| < 1.5 , where W is the •yp centre-of-mass energy, and with fixed-order NLO 
calculations [17] (Fig. 3). The cross sections for both cases are above the predictions, 
in particular for 7] in the forward (proton) direction. NLO resummed predictions [3] are 
closer to the D* data [16], but still too low for high r]. Using different photon parton 
density functions shows some sensitivity to the parton density parametrisation of the 
photon. A tree-level pQCD calculation [18], where the cq state is hadronised rather than 
the c quark, gives a better agreement with the data compared to the NLO calculations. 

D*^ and Associated Dijets 

Events with a reconstructed D*"^ and at least 2 hadron jets ("dijet event") enable one to 
study the photon structure, in particular its charm content. D*^ candidates with p^* > 
3 GeV and \r]^*\ < 1.5 have been selected. The two jets with the highest transverse 
energy, Et, in the pseudorapidity range \ri^'^^\ < 2.4 are required to have E^^^ and 
Ej^ above certain values. The fraction of the photon momentum participating in the 
dijet production, x^^^ , is defined as x^^^ = where y is approximately the 

fraction of incoming electron energy carried by the photon. In LO QCD, direct processes 
have x^^^ = 1 while resolved processes have < 1. Samples enriched with direct 
(resolved) events are separated by a cut x^^^ > 0.75 (< 0.75). 
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FIGURE 4. (a-b) Differential cross sections da/dx)^ ° for dijets with an associated D* in the PHP 

range with Ei^^"^ > 7GeV and Ei^^"^ > 6GeV. In (a), the experimental data (dots) are compared to 
the expectations of the HERWIG MC simulation, normalised to the data, for LO-direct (right hatched), 
LO-resolved (left hatched), LO-resolved without charm excitation (dense hatched) and the sum of LO- 
direct and LO-resolved photon contributions (full histogram). In (b), the data are compared with a 
parton level NLO fixed-order calculation, (c) Ratio of low to high x^^^ for D* dijet events with 
^jeti,2 ^ g 6.5GeV vs. compared to the AROMA and CASCADE MC's. The left-hand point is 
due to PHP events, (d) Differential distributions in \ cosO* \ iox D* dijet events with Ep*^'^ > 5GeV. Data 
(dots) and PYTHIA MC (lines) are shown separately for direct- (open dots/dashed lines) and resolved- 
(black dots/full lines) photon events. All distributions are normalised to the resolved data in the lowest 4 
bins. 



Differential cross sections in x° for ZEUS measurements with C = 37pb^ are 
shown in Fig. 4(a-b) [16] and compared with LO MC simulation and NLO fixed-order 
calculation. The peak at high x^^^ is due to the LO-direct BGF process. The low x^^^ 
tail comes from LO-resolved processes, dominated by photon charm excitation. The 
shape of the x^^^ distribution is in good agreement with the MC simulation with ^ 40% 
resolved contribution (Fig. 4(a)). The fixed-order NLO calculation lies below the data 
at low x'^^^ values (Fig. 4(b)). This could be due to the fact that no explicit charm 
excitation component exists in this calculation. 

The dependence of the virtual-photon structure on its virtuality, Q"^, has been studied 
with ZEUS dijet events containing a D* [19]. The cross section ratio R = a{x^^^ < 
0.75)/cr(x°^^ > 0.75) as a function of is shown in Fig. 4(c) to be approximately 



constant up to ^ 200 GeV^, contrary to the no-charm-tag case, where R falls with 
increasing Q^. The data are compared to two MC models which implement no specific 
partonic structure for the photon, generating all low x^^^ events from parton showers 
in two different schemes. The AROMA [20] model, which implements the DGLAP 
evolution scheme [21], lies below the data. The CASCADE results [14] are much closer 
to the data. 

The angular distributions, dN/d\ cos6'*|, of dijet events containing a D* in the PHP 
regime have been measured by ZEUS [22]. Here 9* is the angle between the jet-jet axis 
and the beam direction in the dijet rest frame. The results are shown in Fig. 4(d) sepa- 
rately for direct and resolved photon events. The resolved processes peak at high | cos 9* \ , 
in agreement with LO MC predictions, as expected for dominant gluon exchange. The 
direct processes are much flatter in | cos 9* \ , consistent with quark exchange. The steep 
rise towards high | cos^*| of the resolved charm events provides evidence that the bulk 
of the resolved contribution is due to charm excitation in the photon. 

EXCITED CHARM MESONS 

P-wave charm mesons (L=l of the cq system) can decay into L=0 states plus a tt or a 
K [23]. They are predicted [24] to appear in two doublets with total angular momentum 
j=3/2 (narrow states) or j=l/2 (broad states). Narrow states, Di(2420) and ^2(2460), 
were observed in the D*tt decay mode and identified as members of the j=3/2 doublet 
with spin-parity — 1+ and 2+, respectively [23]. A charm-strange excited meson, 
Dfi(2536), was found in the D*^K^ final state [23]. 

i:'?(2420) and L>|0(2460) production 

Di and Df* mesons were reconstructed by ZEUS [25] via their decays to D*^nf, 
followed by the D*^ decays, D*^ D^^^s ~^ {K~T:+)TTg{+c.c.). Fig. 5(a) shows the 
"extended" mass difference distribution, M{Kt^t^s'^a) — M{Kt^t^s) + M(D*), where 
M{D*) is the nominal D*^ mass [23] (full dots). A clear excess is seen around the mass 
region of the Z}°(2420) and Z}|°(2460) mesons. No enhancement is seen for wrong 
charge combinations (dashed histogram), where D* and 114 have the same charges. The 
solid curves in Figs. 5(a-b) are an unbinned likelihood fit to two Breit-Wigner shapes 
with masses and widths fixed to the nominal and values [23], convoluted with 
a Gaussian function and multiplied by helicity spectrum functions for = 1+ and 2"*" 
states, respectively. The background shape was parametrised by the form a;" • exp{—/3 ■ 
x + 'f-x'^), where x = M [K t^t^ a) — M{Kt^t^s) — M{ti). The fitted curves describe the 
distribution reasonably well, except for a narrow enhancement near 2.4 GeV (Fig. 5(b)). 
In Fig. 5(c), a similar fit is shown with an additional Gaussian-shaped resonance with 
free mass and width. The fit yielded 211 ±49 entries for the narrow enhancement 
with mass value 2398.1 ± 2.1 (stat.)lo!8 (syst.) MeV. The width was consistent with the 
resolution expected from the tracking detector. The enhancement may indicate a new 
excited charm meson, a result of an interference effect or a statistical fluctuation. The 
number of reconstructed and Df* mesons in the 3-resonance fit are 526 ± 65 and 




FIGURE 5. (a-c) Extended mass difference distribution for D*^ candidates (full dots). The dashed 
histogram is for wrong charge cominations. The curves are the results of unbinned likelihood fits. In 
(a) and (b) the solid curves are a fit to background parametrisation and two Breit-Wigner distributions 
convoluted with a Gaussian function. In (c) an additional Gaussian-shaped resonance near 2.4GeV is 
assumed in the fit. The dotted curves are fitted shapes of the combinatorial background, (d) Effective 
M{D*^K°s) distribution (full dots). The soUd line is a fit to a Gaussian resonance plus background of the 
form A(AM'==^*)'^. 



203 ± 60, respectively. 

The acceptance-corrected fractions of D*^ mesons originating from and 
£'2° in the measured kinematic range were found to be Ruo^w^-k^ /d*± = 

3.40 ± 0.42 (stat.)l°:^^(syst.)% and i?D*o-.D*±^T/D*± = 1.37 ± 0.40 (stat.)+°:^^(syst.)%. 
Extrapolating to the full kinematic phase space by a MC simulation and using the 
partial width ratio, T{Df D+Ti-)/T{Df D*+n-) = 2.3 ± 0.6 [23], the rate 
of c quarks hadronising as D*+ mesons, /(c D*+) = 0.235 ± 0.007 ± 0.007 [9], 
and isospin conservation, the rates of c quarks hadronising as and Df mesons 
are found to be: /(c ^ D?) = 1.46 ± 0.18(stat.)+°:p(syst.) ± 0.06(ext.)% and 
/(c ^ Df) = 2.00 ± 0.58 (stat.)+J:^°(syst.) ± 0.41 (ext.)%. The third errors arise 
from uncertainties in /(c — > D*^) and the Df — > D*^n^ branching ratio. The results 
are consistent with the e+e^ rates measured by CLEO [26]: /(c D^) = 1.8 ± 0.3% 
and f{c^Df) = 1.9 ±0.3%. 

Search for radially excited D*'^ 

Radially excited charm mesons with mass around 2.6 GeV are predicted [27] to decay 
into Dnn or D*nn. A narrow resonance in the D*^n^n^ final state at 2637MeV, 
interpreted as the radially excited D* ^, was reported by DELPHI [28]. No evidence 



for this state has been found by OPAL and CLEO [29]. 

D* ^ candidates were reconstructed by ZEUS [25] from their decays to 0*^^7(^71^. 
No narrow resonance is seen in the extended mass difference M^KimsT^AT^b) ~ 
M{K7nrs)+M{D*). An upper limit of < 2.3% (95% C.L.) is 

obtained in the measured kinematic region for the fraction of D*"^ originating from 
decays within a signal window 2.59 < M{D* ^) < 2.67GeV, which covers theoretical 
predictions [27] and the DELPHI measurement [28]. Extrapolating by a MC simulation 
to the full kinematic phase space and using the known /(c — > D*) value, a D*'^ produc- 
tion limit of /(c D*'+) ■ < 0.7% (95% C.L.) is obtained. A similar 
limit of 0.9% has been reported by OPAL [30]. 

Production of the charm-strange meson Dfi (2536) 

Dfi mesons were reconstructed by ZEUS [5] via the D*'^Kg decay mode with 
Kg — > TT+TT". Kg candidates were identified by using pairs of oppositely charged 
tracks with > 0.2 GeV. A clean Kg t^^t^a signal was extracted after applying 
standard V^-finding cuts. K^g candidates with 0.480 < M{ti2,ti4) < 0.515 GeV were 
kept for the Df^ reconstruction. Fig. 5(d) shows the effective M(D*^Kg) distribu- 
tion in terms of AM^""* + M{D*+)pdg + M{K^)pdg (solid dots), where AM^=^* = 
M{K7r7ig7i37i4) - M{KTnrg)-M (713714) and M{D*+)pdg (M{K^)pdg) is the nomi- 
nal D*^ (K^) mass [23]. A clear signal is seen at the M{Df^) value. The solid curve 
is an unbinned likelihood fit to a Gaussian resonance plus background of the form 
A^AM^""^)^. The fit yielded 62.3 ± 9.3 Df^ mesons. The mass value was found to be 
M{Df^) = 2534.2 ± 0.6 ± 0.5MeV, in rough agreement with the PDG value [23]. The 
last error is due to the uncertainty in M(D*+) pdg- 

The angular distribution of the Dgi signal was studied via the helicity angle, a, 
between the Kg and TTg momenta in the D*^ rest frame. The dN/d cos a distribu- 
tion was fitted to (1 + 7?cos^a). An unbinned likelihood fit yielded R = — 0.53 ± 
0.32(stat.)l°:°|(syst.), consistent with the CLEO value [31] R= -0.23lg;^^. Both mea- 
surements are consistent with R = 0, i.e. = 1+ for the Dgi meson. However, our 
result is not inconsistent with R = —1, i.e. = 1^ or 2+ [32]. 

The fraction of D*^ mesons originating from Df^ in the measured kinematic region 
is Rr)±^D*±K"/D*± ~ 1-77 ± 0.26 (stat.)lQQg(syst.)%. The rate of c quarks hadronising 

si ' 

as D^i mesons, after a MC extrapolation to the full kinematic phase space, is f{c^ 
L>+i) = 1.24±0.18(stat.)lo:o6(syst.)±0.14(br.)%. The third error is due to uncertainties 
in /(c D*^) and the D^^ D*^K^ branching ratio. The rate agrees with the OPAL 
value [33] 1.6 ± 0.4 ± 0.3%. This rate is about twice that expected, assuming /(c 
D\) 2% [25] and 7^ a; 0.3 [8], where 7^ is the strangeness suppression factor in charm 
production. 
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FIGURE 6. (a) Differential cross section in pT^^ for the production of electrons in dijet PHP events 
compared to HERWIG MC shape prediction (solid line). The MC beauty (charm) component is given 
by the forward- (backward-) diagonally hatched histogram, (b) Muon impact parameter, 5, distribution 
for DIS events, with a likelihood fit decomposition into beauty (white region) and charm + fake (shaded 
region), (c) Ratio of measured h production cross sections over theoretical expectation, as a function of 
Q^. The shaded band represents the theoretical uncertainty. 



OPEN BEAUTY PRODUCTION 



The production rates for beauty at HERA are about two orders of magnitude smaller 
than for charm. Theoretical uncertainties are expected to be smaller due to the high 
h quark mass. Enriched h samples have been obtained [34, 35, 36, 37] by studying 
electrons or muons from semileptonic (SL) h decays. Due to the heavy h quark, the 
high tail of the lepton px with respect to the axis of the closest jet, p^?', provides a h 
signature. HI has also exploited, using their micro-vertex detector, the long lifetimes of 
h hadrons to extract h production cross sections by measuring the impact parameter, 5, 
which is the distance of closest approach of SL muons to the primary vertex in the plane 
transverse to the beam axis. 

The ZEUS differential cross section, da/dp^'-, for a PHP dijet event sample ofC = 
38.5pb^^ with identified electrons is compared in Fig. 6(a) [35] with a MC c and b 
production simulation. The shape is fitted by varying the relative contributions of c and 
b. Ab fraction of ~ 15% is obtained, consistent with the MC expectation. 

For a PHP dijet sample with an identified muon, HI has shown [36] that the p'!^'- and 
5 methods provide independent and consistent b production results. The HI analysis 
of a smaller DIS sample from C = lO.Spb^^ [37] uses the combination of the two 
observables. A likelihood fit of bb, cc and fake muon spectra to a two-dimensional 
distribution in S and p'!^'' adjusts the relative weights of all three components in the 
data. It yields a bb fraction of (43 ± 8)%. The 5 projection of this distribution is shown 
in Fig. 6(b) together with the fit decomposition. 

The ratio of the measured visible b cross sections over theoretical expectations [10, 17] 
is shown in Fig. 6(c) as a function of Q'^ [37]. The ratio is roughly constant with and 
the discrepancy between data and NLO calculations is quite significant. The DIS case 
is theoretically cleaner, since at high the resolved contribution is expected to be 
suppressed. 
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